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Randomly Stopped Nonlinear Fractional Birth Processes^ 

Enzo Orsingher^'*, Federico Polito^ 
Abstract 



We present and analyse the nonlinear elassical pure birth process 9i{t), t > 0, and the fractional pure birth 
process 9l'^(t), t > 0, subordinated to various random times, namely the first-passage time Tt of the standard 
Brownian motion B{t), t > 0, the a-stable subordinator s"{t), a G (0,1), and others. For all of them we 
derive the state probability distribution pk{t), k > 1 and, in some cases, we also present the corresponding 
governing differential equation. 

We also highlight interesting interpretations for both the subordinated classical birth process 5\£ (t), t > 0, 
and its fractional counterpart ^'^{t), t > in terms of classical birth processes with random rates evaluated 

Q^ ' on a stretched or squashed time scale. 

p I . Various types of compositions of the fractional pure birth process 7^ " (t) have been examined in the last 

part of the paper. In particular, the processes !Ji^{Tt), 9l^{s°'(t)), ot'^ {T2u{t)) , have been analysed, where 
T2i/(i), t > 0, is a process related to fractional diffusion equations. Also the related process 'H.{s"{T2u{t))) is 
investigated and compared with 91 {T2v{s°' {t))) = 'H." {s°' [t)) . As a byproduct of our analysis, some formulae 
relating Mittag-Lefher functions are obtained. 

Keywords: Fractional nonlinear pure birth processes. Subordination, a-stable subordinator. Fractional 
►^ , derivative. First-passage time, Mittag-Leffler functions, Wright functions, Lamperti law. 

00 

OO ■ 1. Introduction 

(N 

f->^ . We here consider the pure birth process 'H.{t), t > 0, (linear and nonlinear) composed with different 

^^ ' processes like the first-passage time of Brownian motion Tt (possibly iterated n-times) , the sojourn time of 

Brownian motion Tt and bridge &t, and a-stable processes s'^{t). 

The subordination of processes (first introduced by Bochner [J]) has been studied by several authors, over 
the years, in connection, for example, to modelling the wear of instruments during the real working time, or 
is^ " security trading which takes into account fluctuations of the economic activity during the time elapse t (see 

^2 ! Lee and Whitmore Q). 

C^ ' The second part of the paper concerns the subordination of the fractional pure birth process 9{,''{t), t > 0, 

< J^ < 1, with the processes s"{t) and T2a(t), establishes that 5v;''(5"(t)) = 5V;(T2^(5"(i))), and discuss its 
connection with 5V;(5"(T2i,(t))). 

Subordinated processes connected with fractional and higher order partial differential equations are 
treated in numerous recent papers. Most of them concern compositions of time-continuous processes (see 
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for ex amp le Baeumer et al. [21), but also point processes (Laskin |11| . Mainardi and Gorenflo [13], Uchaikin 
et al. [17|, Beghin and Orsingher [3|, Meerschaert et al. 14|). 



Birth processes stopped at different random times can be useful to model branching processes under 
laboratory conditions. For diseases started off artificially, the spread of the infected population can be 
stopped when the experiment leads to convincing conclusions. The cost of the investigation can play a 
certain role in stopping the artificially constructed experiment. The fluctuations of the temperature during 
the effective time t can influence the growth rapidity of cells or of bacteria and thus the population size 
can be thought as a function of the temperature modelled as a random time process. The same reasoning 
underlies experiments in physical studies on chain reactions. In the case of y\C{Tt), where Tt is the sojourn 
time of a Brownian motion on the positive half-line, the experiment can be interrupted immediately (if it 
proves useless), or at the end of the time interval [0,t] (in the case that no evidence can be attained in a 
short time). 

We recall that the distribution of the nonlinear fractional birth process (with one progenitor) reads 

Pr{.^''(t) = k I ^^(0) = 1} = ^^'^' ' ^"=' nf.,,,^„(A,-A™)' ''>'> (11) 

[E,,i{-Xit''), k = l, 

where 

h=0 ^ " 

is the Mittag-Leffler function and Aj,, k> 1, are the birth rates (see Orsingher and Polito 15] ). 
For Xk — X ■ k (fractional linear birth process), formula (jl.ip takes the simple form 

Pr{Af'^(t) = fc| Af^(0) = l}= V f''"M(-l)'"-i^^i(-Amr), fc > 1, t > 0. (1.3) 

^'-^ \m — 1/ 

711=1 ^ / 

For u = 1, we retrieve from (jl.ip and (jl.3p the classical distributions of nonlinear and linear pure birth 
process, by taking into account that Ei^i{x) = e^. 

The simplest subordinator considered is the first-passage time 

Tt^mi{s:Bis)^t}, (1.4) 

where i? is a standard Brownian motion, independent of the birth process considered. For us it is relevant 
that the probability density of (|1.4|) 

q{t,s)ds = PT{Tteds}, (1.5) 



satisfies the following equation 



q2 C) 

^q(t,.s)=2— q(i,s), i>0, S>0, (1.6) 



as a direct check shows. 

In view of (|1.6p we can establish the following relation between the state probabilities 

pUt)= Pi- W'{Tt)^k} (1.7) 

and (fTTj) : 



d^ r°° d 

—p-{t)^~2j^ qit,s)-PT{^''is)^k}ds. 



For v — 1, equation (jl.Sp becomes the second-order difference-differential equation 

^m)^2[\,pUt)-Xk-irk^M^ fc>i- (1-9) 

Furthermore, for v — 1, the probability distribution (11.7^ can be worked out explicitely and becomes 

For < i^ < 1, in light of the well-known integral representation of the Mittag-LefHer function 

i^.,1 -A^ = / 2.^0 . TT^*^' ^^£0,1, 1.11 

TT Jq r^'' + 2r'^ cos t^7r + 1 

we obtain several different representations of the distributions of the subordinated processes. 
For 1/ = 1/2, we have the following result 

p|W = — / ,^2 , ^, PrK^(t)^fcK^(0)^l}d^, (1.12) 

which shows that !K^^'^{Tt) is equivalent in distribution to a fractional pure birth process (denoted by fACwit)) 
with rates Xk ■ W, where V(^ is a folded Cauchy distribution with scale parameter equal to -\/2- 
We have also that 

pl{t)= I Pr{:At(s) = fc}Pr{|C(\/2i)| erfs}. (1.13) 

In other words, 9i^'^{Tt) is also equivalent in distribution to 91 {\C {^/2t)\) , C being a Cauchy process. 
We generalise the previous framework by considering the iterated process 



rpl 



where T} , . . . ^ T" , are independent first-passage times and 



i>0, (1.14) 



r^,+i = inf <^ s: B^{s) = T^+^\ , j = 1, . . . , (n - 1), (1.15) 

where B^ {t) , t > 0, 1 < j < n, are independent Brownian motions. In particular, for i^ = 1 we show that the 
state probabilities 

Pr{^i(0 = fc}=F^(t), (1.16) 

satisfy the 2"th order equations 

^-plit) = 2'"-' {X,~pl{t) - X,-^~pLAt)} ■ (1-17) 

The distribution 91^ (t) =91 (t) (for short) is directly derived and reads 

T-rfc— 1 \ Y^A: e ^ h ^^ ^ 

ii,=iAjZ.™=inf^,,^„(A,-A,„)' ''^'' (1.18) 



For n — > oo, we obtain from (jl.lSp that 

lim;;,(t) = r ' ' (1.19) 

n^oo I 0, fc > 1. 

In the last part of the paper we examine different types of compositions of the fractional pure birth process 
with positively skewed stable processes s"{t), t > 0, < a < 1. For a = v, we show that 

^•'{s-'{t))'^-^{m,), 0<iy<l. (1.20) 

For the stable random variables Si, S2 , it is well-known that the ratio 

>V„^(|)° (1.21) 

(sometimes called Lampcrti law), has probability density equal to 

, , sini^TT r'^^^ 

fwAr)^ 2^-T^ TT' ''>0- (1-22) 

an r a + 2r a cos vk + 1 

i.d. 

Furthermore, we show that 91" {s" {t)) = 'HXT2i,{s'' {t))) ¥" ^ {s" {T2v{t))) ■ We are also able to prove that 

^"(T2„(t)) = ^(r2,(T2„(t))) =■ ^(r2,„(f)) = ^""(i). (1.23) 

As a byproduct of our analysis we obtain the following integral relation between Mittag-Leffler functions of 
different indices: 

sinz^TT f°° r"-^ / i \ 

E,o.A-Xrntn = / 2. , o . —E^,J-rX^e)dr (1.24) 

TT Jo ^ + ^"f" cos t/TT + 1 V / 

sinaTT f°° r°'^^ / 1 , \ , 

= / 2a ^o a --E,,J-rX^ntndr, 0<a,u<l. 

TT Jo ^ + 2?" COS aTT + 1 V / 

2. Subordinated nonlinear birth processes 

In this section we study in detail the nonlinear pure birth process stopped at Tt and we derive the state 
probabilities p — Pr{.^(r4) = fc | 5v;(0) = 1}, fc > 1, and the corresponding governing differential equations. 

We give some information about the process Jt{t), t > 0, evaluate explicitly its mean value E5v;(i), and 
discuss also the linear birth process (sometimes referred as Yule-Furry process). 

2.1. Preliminaries 

The state probabilities pfc(^) = Pr{5v; (i) = fc | 5v;(0) = 1} read (see e.g. Gikhman and Skorokhod 7], page 
322) 

[e-^i*, fc = 1, i> 0. 

For the case of uq progenitors (see Chiang [Sl, page 51), formula (|2.ip must be replaced by 

rfc-l 

p^{t) = <; ^^^="'' ''^■^™=»o nt„„,^,„(A,-A„)' "' - (2.2) 



fc = TT-O; i > 0. 



-A„ot 



We assume that ^^. 1/Afc = oo in such a way that the process is non-exploding (see Feher [6|, page 452). For 
a discussion on this point, consult Grimmett and Stirzaker [8[, page 252. The probabilities (|2.ip satisfy the 
following difference-differential equations: 

-jTPkit) ^ -\kpk{t) + Xk-ipk-iit), k>l. (2.3) 

We have our first result in the next theorem. 
Theorem 2.1. The mean value of the nonlinear birth process is 

oo I k k . I 

fc=l [ m=ll=l,ly^m "* j 

Proof. From equation (|2.3p . we have that 

oo , oo oo oo 

E ^i;p'=(*) ^ ^ E ^^kpk{t) + E k\k-ipk-i = E ■^kPk{t)- (2.5) 

fc=l fc=l fe=2 /i:=l 



By integrating both members in (0,t), we obtain 

oo OO' -.i 

E fcpfc(i) - 1 = E '^'^' / Pk{s)ds 
Jo 



fc=l fe=l 



t oo I „t fc-1 fc _\ 

— I / — — g A„ 






'^■-1 fe ^1 -A,„A 



l-e-^* + E|nA.E ^'"^' 



oo k k -. 

i--^^*+EE n 3^(1--^^*) 

fc— 2 m— 1 /^l,//m 



fc=2 I rn=ll=l,l^m 



I p-A,„t 



OO I A: fc 

E i-E n -^e 

fc=l I m=ll=l,l^m 



A; — A„ 



(2.6) 



so that_formula (|2.4I) emerges. In the second-to-last step of (|2.6p . we applied formula (3.12) of Orsingher 
et al. jl6[ and, in the last step, we considered that, for k = 1, the set of numbers {1 <l <\, I ^ m ^ 1\, is 
empty and the 

n jtk.- '"> 

is taken equal to 1 by convention. D 



Remark 2.1. As a check we can extract, from ()2.4|) . the mean value in the linear case Am — m ■ X. Since 



k k . , k 



^^ „-Amt_Y- l...(m-l)(m + l)...fc 



V TT ^-Amt ^ Y^ 1... V/>t- iyv»t-M;...fi. ^^„t 

^ -l-J- AZ-Am ^ fm-l)(m-2)...l-(-l)"-i-l...(fc-TO) ^ '^ 



m=l l^l^l^m 



^ (m - l)(m - 2) ... 1 • (-1)"-! • 1 . . . (fc - to) 



m=l ^ ' m=l ^ ' 

we have that 

1-E n A[^^"'"*=(i-^"*)'- (2.9) 

m=l (=l,;:^Tn 

From this we readily have that 

EN{t) = 1 + ^(1 - e-^')" = 1 + ^_ , - 1 = e^*. (2.10) 

The aim of this section is to compose the process 5v;(i) with the first-passage time Tt — inf(s: B{s) = t), 
where i? is a Brownian motion independent of ^ (t) . 

Remark 2.2. The probability density of Tt — inf {s: B{s) ~t},t> 0, where B(t) is a standard Brownian 

motion, namely 

_£. 

Pr{Tteds}/ds^q{t,s)^t^^, (2.11) 

v27rs-' 

is the solution to the Cauchy problem 

U(0,s) = <5(s) 



9^q{t,s)=2f^q{t,s), t>0,s>0, ^2^2^ 



as a simple check shows. 

2.2. Pure birth process stopped at Tt 

Theorem 2.2. Let ^{t), t > be a classical nonlinear pure birth process and let q[t,s), s > 0, i > 0, the 

law of Tt. The process 9i{t) = •^c{Tt), t > 0, has the following distribution 

hit) = I ^^'=' ^' ^"=' mJ,^^(>^7-K.) ' ^>l'^>0' (2.13) 

g-tv2Ar^ fc = 1, i > 0, 



ana mean 



value equal to 



A; 



E^(r,)=i+5: 1-^ n Y^'""^]- (2.14) 

fe=l \ rn=ll=l,l=i,n^^ '^"^ 



The distribution (12. 13^ is non-exploding under the condition that ^j, l/Xk — oo. 



Proof. The state probabilities are derived by straight calculations and by resorting to the Laplace transform 
of q{t, s) which reads 

e^^'q{t,s)ds= r e-"^^^ds = e-'^^. (2.15) 

Jo v27rs3 

We treat the case fc > 1 as follows. The case /c = 1 is analogous. 

/•OO 

Pk{t)=Pr{^{Tt) = k\^{0)^l}= pk{s)q{t,s)ds (2.16) 



A; — 1 k —Xs 4- —— ^'~1 ^ 1 

1~T \ \~^ ^ "^ '^^ J 1~T \ \~^ —t\/2\ 

In view of Theorem 12.11 we can evaluate the mean value 

^oo , _ii OO / fe fc , \ 

u 

In the linear case (j2.14p can be written as 

OO fc 

E5V:(Tt) =^^(-l)™e~*'''^A^. (2.18) 

fc=0 m=0 

On the other side, this sum diverges because 



E^ 



(Tt) = / e^^PrjTt e ds} = cx). (2.19) 

Jo 



Remark 2.3. Note that Vi, pk{t), k > 1 is a proper probability distribution because of the composition 
"Kit) — 9(XTt). The process can be appropriately interpreted by rewriting (|2.16p as follows 

fc-i k „oo y- - — 

h{t) = X{XjY.T^k 7^-^-J ^"'""^^^ (2.20) 

-->fc-l fe -, _J_ 

The process 'K{t), t > can be viewed as a classical nonlinear pure birth process evaluated at time t^ with 
random birth rates QXk, k>l, where Q is an inverse Gaussian random variable with p.d.f. 

The composition oi 9i{t), t > 0, with Tt leads to a second-order time derivative in the governing equations, 
as shown in the next theorem. 



Theorem 2.3. Let pk{t), t > Q, k > 1, he the distribution of the process ^{t) — 9l{Tt), t > 0, where Tt 
is the first-passage time process of the standard Brownian motion, having transition density q{t,s), s > 0, 
t > 0. The state probabilities pfc(t), t > 0, fc > 1, satisfy the following difference- differential equations 

^Pfe(i) = 2 [Xkhit] - Xk-ih^iit)] , fc > 1, (2.22) 

where Xk, fc > 1 are the birth rates of the nonlinear classical birth process 'H.{t), i > 0. 
Proof. Since 

POO 

h(t)= / pk{s)q{t,s)ds, (2.23) 

Jo 

by taking the second-order derivative w.r.t. t, in view of Remark 12. 2| we have that 

^2 f°° d = f°° d 

— Pfc(t) = 2y pk{s)-^^q{t,s)ds^2q{t,s)pkis)\l^^-2j —pk{s)q(t, s)ds (2.24) 

/•OO 

= -2 / q{t, s) [~XkPk{s) + Xk-ipk-i{s)] ds ^2 [Xkpk{t) - Xk-ipk-i{t)] ■ 
Jo 

In ()2.24p . we considered that pk{0) = 0, for fc > 1. D 

Remark 2.4. In the linear case, some calculations suffice to show that 

" /fc-1 

m=l ^ 

and the state probabilities satisfy the equation 

d^ 

-^Pk{t)=2Xpk{t)-2X{k^l)pk-i{t), fc>l. (2.26) 

2.2.1. Iterated compositions 

Theorem 2.4. Let !At(i), t > 0, be a classical nonlinear birth process. Let T^ , T^,...,Tp, n G N, 6e 

first-passage times of n independent standard Brownian motions. The process 



Pk{t) = Pr{N{Tt)^k}^y^{ (_i)™-ie-*V^^, fc > 1, t > 0, (2.25) 

^ — ' \ TO — 1 / 



has the following distribution 



lrp2 



t > 0, (2.27) 



^k^i . ^k „-t^i'2i' ^) 



, miiA,El^. £^;^(.,_.„) , fc>i,i>o, 



Ffc(i) = <(^^^=' '^"=' nr=i,M™(^'-^'") (2.28) 

e~^^r2i'-^\ fc = l, i>0. 



Proof. We start by proving the case n = 2 since the case n = 1 is aheady proved in Theorem 12.21 We omit 
the details for the case fc = 1 and directly treat the case fc > 2. We have that 

Pk{s)q{t,s)ds^ l[X,J2-^—— —^==ds (2.29) 



fc-1 



= 1 Ul^lA^m i^l ^ ^rn) Jo 



fc-1 fc 



-,-sV2A,„ __^^ 



V2Trs^ 



1 



-t\Jl^f2\^„ 



j = l m=l ll/=l,;/m V'^' '^'") 

It is now straightforward to generalise formula (j2.29p for n compositions, as follows 

fc-1 fc 



Pfc 



w^n^.E 



1 



1 ^n 1 



'tA 2" 2 



rfc 



fc-1 fc 



(2.30) 



n^^E^ 

j = l m=l ll/=l,;5^m l'^' '^m) 



-tAr"2('"2^) 



When n ^- oo, equation p.28p becomes 
lim pfc(i) 



e-"n^i^A,E:^i nuJ(A,-A„o =Q' ^>i' 



-2t 



fc = 1, 



n 



(2.31) 



because of formula (3.4), page 51 of Chiang [5|. Therefore, the process (j2.27p can either assume the state 
k — 1 with probability e~^*, or explode with probability 1 — e~^*. 

Theorem 2.5. Let pk{t), t > 0, k > 1, be the distribution of the process 



■K{t) = 'H, 



Irp 



rp2 



t >0. 



The state probabilities pk{t), t > 0, k > 1, satisfy the following difference- differential equations 

-rp^Pkit) = 2^ ^ {Afepfc(i) - Afc-ipfc-i(t)} , 
where Xk, k > I, are the birth rates of the nonlinear classical birth process "Kit), t > 0. 
Proof. For n = 1, equations (|2.33p reduce to equations (|2.22p . For n = 2 we have that 



(2.32) 



(2.33) 



di4 



~Pk{t) - 



OO t^OO 



'pk{wi)q{w2,wi)j—q{t,W2)dwi dw2 



(2.34) 



"'0 

2^/ / h{wi)-p-^q{w2,Wi)q{t,W2)dwidw2 
Jo Jo Cw;2 

(•CO coo a2 



Pfc 



JO 

OO /"OO 



d 

{wi)^-^q{w2,wi)q{t, W2)dwi dw2 
OW2 

= 2''/ / Ffc(wi)T; — q{w2,Wi)q{t,W2)dwidw2 
Jo Jo owi 

/•OO /"OO 7 

= -2^ / 3 — pk{wi)q{w2,Wi)q{t,W2)dwidw2 

Jo Jo dwi 

= 2 |Afepfc(i) — Afc_ipfe_i(i)| . 
The above reasoning can be generalised, thus arriving at equation p.33p . 



D 



2.3. Other compositions 

In this part we present the distributions of the classical nonlinear birth process 5\£(t), i > 0, stopped at 
various random time processes, namely the sojourn time T^ of a standard Brownian motion, the sojourn time 
©( of a standard Brownian bridge and the stable subordinator 5"(t) of order a G (0, 1]. 

We start first by considering the nonlinear birth process at time 

Tt = / I[o,oo){B{s))ds = meas{s < t: B{s) > 0} . (2.35) 

Jo 

The process !AC{Tt), is a slowed down birth process. In the next theorem we provide its distribution. 
Theorem 2.6. We have that 

fn^^^A V'^ '^'^"'"^"(^^"O fc>lt>0 
Pr{at(rt) = fc}= <^^J=i^^™=inf=i.M™(^'-^'")' ''^^'^^"' (2.36) 

[e-t^i/o(|Ai), k^l,t>0, 

where 

oo 2/i 1 

k=Q *■ ^ 

is the zero- order Bessel Junction with imaginary argument. 

Proof. The derivation of (j2.36p is based on the evaluation of the following integral: 



L 



-s\,„ '^S .-iA ^ /"t 



-sA,„ ^ e^^A™/ x„ . (2.38) 



Remark 2.5. In view of the integral representation of the Bessel function 

Io{z)^^ f\^^°^^dd, 
we can give the following alternative, interesting representation of (|2.36l) . 

Pr{^{Tt) = k] = ^ j Prl^i (isin2 =k\dd. 



In other words 

where & is a random variable uniform in [0, 27r] . 
Theorem 2.7. For the nonlinear birth process stopped at 



'0 

B{s), s > 0, being a Brownian bridge, we have that 



Alt 



fc = 1. 



D 



(2.39) 



(2.40) 



^{Tt)'^-^(tsm'^\, (2.41) 



®t= / /[o,oo)(-B(s))ds, (2.42) 



Pr{^{&t) = fc} - I ^'^* 1\ ^™=i ntiJ#m (a,-L) ''' ^'"7 ' fc > 1' (2.43) 



10 



Proof. The calculation 



fe-i 



Pr{^(0O = fc} = n ^. E TT. 7—^- 

1=1 [U=lMrn {■^l ^ -^"1^ 



j = l m= 



I e^^-'Prjet eds}, (2.44) 



is sufficient to arrive at result (j2.43p . once the well-known fact that (I2.42p is uniformly distributed in [0,t] is 
considered. D 



Remark 2.6. For the linear birth process, the distribution (|2.43p takes a very simple form as the calculations 
below show. Since for Xk = X • k, X > 0, we have that 



n 



A, 



l=lj^m 



Xi — A^ 



{~ir-\ 



(2.45) 



we can write that 



Pri 



{N{e,) = k} = ^(i-Y.(t) i-^Y"-'^-""'' 



(2.46) 



T, t f!)<-') 



Xkt ' — ' \m 

m=0 



\7n ~\mt 



l-e 



-\t\ 



Xkt 



k>l. 



The distribution (J2.46I) is logarithmic with parameter 1 — e ^* . In the logarithmic distribution with parameter 
< q < 1, we have that 

EL = 



(l-g)log(l-g)' 



YarL 



(1 - <Z)2 log(l - q) 



log (1 - q) 



In our case q — 1 — e so that 



^\t 



EiV(e5t) 



At 



WarN{<&t 



e^*(e^* - 1) 



Xt 



l-e 



-At 



Xt 



For large values of t we have that 



,At 



WarN{(3t) 



EN{t) 



Xt Xt ' 

gAt(-gAt _ i^ YarN{t) 



Xt 



Xt 



(2.47) 
(2.48) 

(2.49) 
(2.50) 

(2.51) 

(2.52) 



Theorem 2.8. For the nonlinear birth process stopped at an a-stable time S°'{t) with distribution qa{t,s) 
and Laplace transform L e~^''qa{t,s)ds = e~^^ , we have that 



Pr{:^£(5"(^))-fc} 



lli = l Aj l^rn=l fP 
■^ 1 ii = l, 



,(A,-A„)' 



e-*^?. 



k>l, 
k = l. 



(2.53) 



11 



Proof. The following calculation is sufBcient to prove result ()2.53p : 

fe-l k 



Pr{^(5"(t)) = fc} = n ^^- E TTk 7\ r^ / e-^"^g„(i,s)ds. (2.54) 



J — 1 m- 



n 



Remark 2.7. Formula (|2.53p can 6e further worked out as follows. 

PrMs-it)) ^k} = Ha, E n. '"'"'_, , (2-55) 

by exploiting the self- similarity of s"{t) 

K 1 k -. pOO 

= n ^^ E TTk 77 T^ / e-^-H-^q4l,t-^s)ds 

j = l m=l [h = l,l^my'^l ^ -^rn) -Jo 

k — 1 k -. ^oo I 

= n ^. E 7^^ ^— T / ^-'-''''lail^OdC. 

j=l m=l i.h = l, l^rny/^l- ^ -^rn) JO 

The last result implies the following representation: 

^{S"{t))=-9{,{t^z), (2.56) 

where Z has distribution ^^(IjC), (^ > 0. 

Remark 2.8. // we assume a = 1/2" in the first line of (|2.55p . and s = t2^^2^ in (|2.28p . the distribution 
(|2.55p .suggests the following unexpected relation: 

Pr^9i(s^{t)) =k\ = Prl9l T^2 =k\, fc > 1. (2.57) 

Remark 2.9. Many other compositions can be envisaged and in some cases they provide curious results. For 
example, we consider the standard Cauchy process C{t), with law h{t,s), i > 0, s G M, satisfying the Laplace 
equation 

We can show that 5v;(|C(t)|), t > 0, is a birth process whose state probabilities p1{t), t > Q, satisfy the 
difference- differential equations 

^Pm = ->^lpm + Afe-i (A, + Xk-i)pU{t) - A.-iA,_2P,t„2 W- (2.59) 

3. Subordinated fractional birth processes 

In a previous work of us (see Orsingher and Polito [15|) we constructed and analysed a fractional (possibly 
nonlinear) pure birth process 5\£'^(i), t > 0, v £ (Oi 1] by exchanging the integer-order time derivative with 
the Dzhrbashyan-Caputo fractional derivative in the difference-differential equation (j2.3p governing the state 
probabilities. We recall that the Dzhrbashyan-Caputo derivative has the form, for < i^ < 1 

^m) = /nT^^o(^'^^' '^<^<^' (3 1) 

dt--'" [fit), v^l. ^■' 
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In this section we examine properties of the subordinated processes 5V;'^(Tt), i > 0, 9i'^ {T2p{t)) , smd 9i" [s"' {t)) , 
i > 0, i/, a, /3 e (0, 1], bringing to the fore some interesting relations and discussing the interpretation for the 
resuhs obtained. 

3. 1 . Preliminaries 

The state probabihties v'^{t) = Pr{5v;'^(i) = fc}, fc > 1 of the fractional pure birth process have the 
following form 

rf«)4n-'''^''-nftS«^' '^ '■'""■ (3.2) 

i£^^,i(-Ait''), A; = l, t>0, 



where E,^^i{—(t'') is the Mittag-Lefher function defined as 



Eu,i{~Cn^Y. v( hlw CeK, ^>0, (3.3) 



h=0 

and with Laplace transform 



'E,^,i-Cndt=^—-, iy>0. (3.4) 

z + (, 



A useful integral representation for Eu^i{—(^t^) reads 



2r^ cos z/TT + 1 



E.A^n = —^ I 2. , o . dr, V e (0, 1). (3.5) 



In a previous work (see Orsingher and Polito [15|) we proved a useful subordination representation for the 
fractional pure birth process p.2p . This can be viewed as a classical birth process stopped at a random time 
Tivify possessing density function coinciding with the folded solution to the fractional diffusion equation 

g(x,0)=<5(a:), 

with the additional condition gt {x, 0) = for 1/2 < i/ < 1. In other words 'K.'^{t) — 9i {T2uii)\ i > 0. It can 
be shown that fT2u{s,t) = Pr{T2i,(i) G ds} is also a solution to 

d^f df 



(see Orsingher et al. |16l|). 

Theorem 3.1. The fractional nonlinear pure birth process is a renewal process with intermediate waiting 

Pr{T^ e ds} = Xks''-^E^^^{-Xks'')ds, fc > 1, s > 0, (3.8) 



times Tj^ with law 



where T^ is the random time separating the kth and (fc + l)th birth. 

Proof. We prove this result by induction. Denoting Z^ = T^' + • • • + T^, we can certainly write that 

Pr {rr + • • • + T;: edt}= f Pr {Ti: e d(< - s)} Pr {T^ + • • • + T^i G ds} , (3.9) 
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where Pr {Tf + • • • + Tj^_i G ds] /ds = ^Pr {'H.'^{s) > k}. By resorting to Laplace transforms, from p.9p . 
we obtain that 



/ e-^*Pr{Tr 
Jo 



T^^ edt}^ / e-^*di / Pr{T^r e d{t - s)}Pr{Ti'' + • • • + T^^i e ds} (3.10) 

/•OO /"OO 

/•CXD /'OO 

Jo Jo 



n/ e--Pr{i;''eds} = n7I7 



A. 



+ A, 



We observe that 



and that 



PrlTi" e ds}/ds = ^Pri^-'is) > 2} = -^E^.^i-X^s") = Ais''-iS,,,(-Ais'^), 
as as 



Pr{Tr + Ta"' G ds} = -^ [1 - Pr{^''(s) = 1} - PrK^(s) = 2}] 
as 



(3.11) 



(3.12) 



Ais'^-^£;,,,(-Ais'') + Ai 



' S,,,(-Ais'') , ^ £;...(-A2s'') 
M — ^ ^ \- M 



A2 — Ai 



Ai — A2 



A2 ^ Ai 
For k = 2, relation p.lOp simplifies to 



^'^' s"-' [E,,u{-Xisn - E,.A~^2sn] 



/ e"^*Pr {T'( + T^ e dt} 
Jo 



Ai A2 



fi'^ + Ai fi" + A2 
and this coincides with the Laplace transform of p.l2p . 
Theorem 3.2. The mean value ¥,9l'^{t), for the fractional nonlinear pure birth process has the form: 

k k 



E^-^w^i+E 1-E n 



A, 



A; — A„i 



--E',y.l(— Amt") 



k—l I 771 — 1 t — l^l^m 

Proof. In light of the subordination relation fH.'^{t) = 'K{T2v{t)), and of result (|2.4[) . we can write that 



¥.9i^t)^ E^(s)Pr{r2,(i) eds} 

Jo 



(3.13) 
D 

(3.14) 

t 
(3.15) 



00 



k k 



k—l I m—ll—l,l^7n 

00 I k k 

i+E 1 E n 

k=l I m=ll=lj=£m 



A, 



A/ — At,! Jo 

A; 
A/ — At 



/>CXD 

/ e-^-^Pr{T2,(i)eds} 

Eu^l{—Xmt'^) 



n 
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In the previous steps we assumed that 

/ e-^-'fT,At,s)ds= e-^^-'Pr{T2At)eds} = E,s{-\mn- (3.16) 

Jo Jo 

We give here some details of this resuh. The density fT2„{z,s), z > 0, s > 0, is obtained by folding the 
solution of the fractional diffusion equation 



d^-'u (92 



which reads 



(3.17) 



z" \ z 

where W-u,i-u{—^) is a Wright function defined as 



u{z, s) = — VF_,,i_, (-4) ' s > 0, z > 0, (3.18) 



r=0 

Therefore 



^-^^ r\i (1 — v(r + 1)) 






oo 



_ 1 f> (-1)" 1 _ 1 ^ _ ( i_ 

z" ^ X'^^Ti-i^k + 1 - J/) (z")'= A^z'' "'' " ^ A^z'' 
(by formula (5.1) page 1825, Beghin and Orsingher [3[) 

Remark 3.1. We can extract, from p.l4p . the mean value of the fractional linear birth process obtained 
in Orsingher and Polito ]la] . formula (3.42), as follows. By considering that X„i = X ■ m, formula p.l4p 
becomes 

°° f '= i,\ ) 

k=l { m=\ '^ '' ) 

= i + E|i + E(-ir(!^)^M(-Amn| 

oo h / , \ OC k y J \ 

= i + EE L (-ir^M(-Amn = EE L (-ir^M(-Amt'')- 

k=lm=0 ^ ^ fe=Om=0 ^ ^ 

/ti order to obtain the desired result we pass to Laplace transforms and extract from p.2ip that 

/ e-''*EiV''(i)dt = E E ("!)''/ e-^'*E^i{-Xme)dt (3.22) 

fc=Om=0 ^ ^ ' •'" fe=Om=0 ^ ^ 
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oo 



fOO p-Wfl" 



Jo M'' - A 

By inverting the Laplace transform above, we can conclude that 

EN''{t)^E^,i{Xt''), (3.23) 

thus confirming our previous result. 

Here we remark that another interpretation in terms of random birth rates can be highhghted. If we write 



•^{t) = Pvi^'^it) = k I 9l-{{)) = 1} = r Pk{s)PT{T2,{t) e ds} (3.24) 

Jo 

K Ik ^ p QQ 

II^.Et^i 7^-^-J e-'-'P^-{T,At)eds} 



3 = 

fc-1 fe 



1 r°° . ..I... / s 






fe-1 fc 



n ^. E T^i ^— T / e-"«*'v_M-.(-Orfe 



oo 



W^-^.i-.l-OPrl^lt") = k I ^(0) = l}d^, 



we have that a fractional nonhnear pure birth process can be considered as a classical nonlinear pure birth 
process evaluated at a rescaled time f^ and with random rates AfeS, fc > 1, where S is a random variable 
with density function 

/H(0 = W^-.a-.(~0, eeK+. (3.25) 

From (I3.24I1 . the following interpretation also holds: 

^■'it)'^-^{Ef'). (3.26) 

Note also that, from p.Sp and p. 241) . we have that 

E.A-^n= e— * W^_.,i_.(-r)dr - / e^ * ^ —dr, (3.27) 

Jo T^ Jo r''^ + 2r'^ COS VTT + 1 

which illustrates an interesting relation between the Wright function and the law of Wi (see (|1.2ip ). Equation 
(|3.27p can be derived directly as follows. 

Jo At- "^''"^V At-y At- £-m!r -^ + 1-^ (At-)"7o 



rn— 



1 ^ (-1)™ 1 _ 1 f 1 



A*'' ,^0 r (-^ + 1 - ^) 7 h^t'^r l^t'^ ' V 7At- 

(by formula (5.1), page 1825, Beghin and Orsingher [3|) 
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E^^ii-jXt") 



g-(7A)i'tr^iy-l 



"■ Jo [r'^ + COS utt) + sin un 



-dr. 



This yields ([3:27| for z = 7A. 

For further details on fractional pure birth process the reader can refer to Orsingher and Polito ^ISj while 
Mittag-Lefher functions are extensively analysed in Kilbas et al. |9j. 



3.2. Fractional pure birth process stopped at Tt 

We consider here the composition of a fractional nonlinear pure birth process, denoted as 9iC'^{t), t > 0, 
V e (0, 1] with the first-passage time Tt of a standard Brownian motion. In the following theorem we derive 
an interesting integral representation for the state probabilities p^(t), t > 0, fc > 1, oi^t'^it) = 9i'^{Tt), t > 0, 

i^e(0,i). 

Theorem 3.3. Let 9i'^(t), t > 0, z^ G (0; 1); be a fractional nonlinear pure birth process and Tt be the 
first-passage time process of the standard Brownian motion with distribution q(t,s). The state probabilities 
Pki^) ~ P'"' {'K'^ {Tt) = k I 5v;''(0) = 1} possess the following integral form 



rfc-l 



llj = l ^i l^m=l Yl''_ 



1 1 roo -E2i-,i(-a;^''e""^)--E2i-,i(-a;^''e ""'') 



x+ty 2\^ 



dx, fc > 1, 
fc= 1. 



Pkit) = S 1 ,00 g2.,l(-:r^-e--)-i;2.,l(-x^-e— ) , " , , (3-29) 

I r 

x+ty2\i'' 
Proof. It is sufficient to prove p.29p in the case fc > 1, since the case /c = 1 is analogous. We have 

fc-l k 



Ffcl 



kit)- / p]:{s)qit,s)ds = '[[X,Y. 



1 



j = l m=l Ill=l,l^m i^l ^ ^m) "'0 



te 2s 

E,y^i{-X„is'') - ids. 

V2^ 



(3.30) 



In order to prove p.29p . by taking into consideration formula p. 51) . we do the following calculations 



E^^ii-Xms") . -ds 



\/2^s3 



te 2s 



smi^TT r e 



f—l^—rX^s 



V2^ Jq tt r2'^ + 2r^ cos i^n + 1 



dr ds 



(3.31) 



smz^TT 



TT 7o ^^^ + 2r'' cos VTT + 1 Jq 



- te^ — 



V27rs3 



■ds dr 



^v-l 



ty2r\^ 



y.2V l_ 0„y 



2i7r 



1 



2r^ cos j/TT + 1 
1 



dr 



^iy-lg-tV2rA^. 



I^ I /^iTTI/ ,»,I^ 



r" + e 
1 



^i^-lg-tV2rA^ 



ITT ./o V w^" + e"" w^'' + e""" 



^2iy-lg-t«;V2A^ 



TT Jo (r'' + e*'''')(r'' + e-*'''') 



-dr 



-twV 2Xi;, 



dw. 



By using the Laplace transform p.4p we obtain 



-E,., 1 ( - A,„ s "^ ) ^== ds 
v27rs3 



(3.32) 
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1 



e-"'^^2^,i(-x^''e"^)da;- / e-'"''E2^A{-x^''e-''"')dx 



-twV 2\f, 



dw 










da;. 



x + t\/2Xik 
Formula p.29p is then proved by combining p. 301) and (I3.32p . 



D 



Remark 3.2. If Xk = kX, k > I, the state probabilities p^(i) = PrlN^iTt) = k \ 7^(0) = 1} of a fractional 
linear pure birth process stopped at Tt read 



m=l Vm-1/ V -L; ijr Jo 



tV2\¥\ 



n'^(i\ — I a;+tV2Ai'r?x 

PfcW " S j_ .00 g2..l(-^"-e""-)-i;2.,i(-:r^-e-"--) 

ITT Jo ' ~ '^•^' 



2;+tV2A>^ 

T/iis result can be obtained by means of methods similar to those of Theorem 



dx, fc > 1, 
fc= 1. 



(3.33) 



Remark 3.3. By considering formula p.30p and the representation p.Sp we can give an interesting inter- 
pretation of the process H.^''^{t) = 9i^''^{Tt), t > Q, as follows (again, we treat the case k>l since the case 
k = 1 is analogous) 



Vk{t) = / Pkis)q{t,s)ds 
Jo 

fc-1 fe 



(3.34) 



7=1 m=l [[l=l.l^m i^l ^™) -^0 






ZG 2s 

E^^i{-X,ns'') . - ds 



Li ^ 21^ T-rfe rx? - A "l TT Jn r'^'^ + 2r'^ cosvTT + 1 



ty 2r\^ 



dr 



k-l k 



j=l m=l lli=l, 



t\j2r\^ 



"■ 7o ^^"^ + ^J''^ COS V-K 

If V = 1/2 we obtain the following expression 

. -i/'OO -1 k—1 k 

V2 



= l,i/m 



(A; — A„i) 



dr. 



Ffc 



,i#m (V2f A, - V2f A™) 



-t\/2rA„ 



dr (3.35) 



1 



-Pr{5v;„(t) = fc| 5v;^(0) = l}dw 



where 9iw{i), t > 0, is a classical nonlinear birth process (|2.ip with random birth rates (WXk), fc > 1 where 
W is a folded Cauchy r.v. with p.d.f 



fwiw) 



V2 



-(^ + 1) 



w e 



It is possible to highlight a further interpretation by rewriting formula (|3.35p in the following way 

k-l k 



P'kit) 



JTy f'OO 1 '^^ J- _ _ 1 



(A/ — Am) 



3-"^'"ds 



(3.36) 



(3.37) 
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= f Pr{9i{s) = fc}Fr| C{V2t) e ds\ 
where C{\2t), t > is a Cauchy process with rescaled time, possessing transition density 

i > 0, s € R. 



fc{t,s) = -- 



^.s2 + (V2i) 
The process •K^^'^it) = 5v;^/^(Tf) can thus he written as 9l'^/'^{t) = •K{\C{\/2t)\). 



(3.38) 



3.2.1. Iterated compositions 

In the next theorem we present the explicit form of the state probabihties p'jl{t), t > 0, k > 1, for the 
process 



n.''{t)^^-' 



rr.1 



rp2 



t>o, ve (0,1] 



(3.39) 



and in the following remark an interesting interpretation for that process when v — 1/2", rt G N, is given. 

Theorem 3.4. Let 9i'^{t), t > 0, v E (0,1], be a fractional nonlinear pure birth process and let T^ , T^ , 
...T", n e N, be n independent first-passage time processes at t of the standard Brownian motion. The 
process 



^•'{t) = ^'' 



rpl 

lrp2 



has the following distribution 



Kit) 



nfc— 1 \ Y->te 1 ; 

1 = 1 ^j Z^m=l fF (\,-\ ~l ' 

Sin i^TT r -^^ re i ^ 

71 Jo ^2(y I O..M „„„ .._ I 1 n 



111/ 7r rO' 

71 Jo 



t>0, ly e (0,1] 



r e ^ 

r-^'^+2r'^ cos u-n+l 



(3.40) 



'+2r'^ cos 1/71+1 ' 



dr, fc > 1, i > 0, 



fc= 1, i>0. 



(3.41) 



Proof. We start by proving the case n — 2 since the case n = 1 is already proved in Theorem 13.31 We omit 
the details for the case fc = 1 and directly treat the case k > 2. We have 



Pkis)q{t,s)ds 



(3.42) 



fc-i 



n^.E 



1 nti.;#m (^/ - ^m) "" ^0 r2- + 2r-cosz/7r + Uo 



„v-i 



„v-l 



j = l m- 
fc-1 fe 

ii J Z^ -prfc (- \, _ \ ^ TT /n r^'' + 2r'' COS z/TT + 1 ■ 



-.(2.)iA#i£lfl^^^^ 



>/27rs"" 



t'\/2V2rA^ 



dr. 



It is now straightforward to generalise formula p.42p for n compositions, as follows 

fe-l k 



Pk 



m-U^sE 



1 



sm VTT 



fc-i 



1 m=l I\l=ll^m i'^l - -^™) '^ "'O 



f.2i/ I Ot^z/ 



2r'' cos vn + 1 



smi/TT 



J Z^ -prfc CAi - A ") 71" Jn r'^" + 2r'' cos j/tt + 1 



g-tr^A^2^?-i2.^^ (3.43) 



p-tr2^A#"2('-5V)^^^ 



D 
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Remark 3.4. Analogously to Remark \3.3\ for v = 1/2", n G N, zt is possible to interpret formula (I3.41|) as 
follows 



fe-i fe 



TT /"CX) 



1 m=i Ui=i,i=i,n (^i - ^m) '^ -'o r^^^ + 2r^ cos ^ + 1 



1 _1 



-tr 2^2 



(-*); 



•dr (3.44) 



(r = 2/^") 



7r />oo 



dr 



^ Jo ^^ + 27" cos 2 
Therefore, the following representation holds 



■^^-^Pr{^ {tr2^-^) -fc} 



5^2" 

where fl is a random variable with density 

h{r) 



1 rp2 



^ (tr!2i-*), 



^ r2 + 2r cos ^ + 1 ' 



re 



The density is a unimodal law which, for n — >■ 00, becomes 

fir) ^ 



(l + ry 



r G 



(3.45) 



(3.46) 



(3.47) 



3.3. Fractional pure birth process stopped at S°'{t) 

We consider the fractional nonlinear pure birth process stopped at a stable time 5"(i) of order < a < 1 
with Laplace transform 

Ee-«°(*) = / e^^'*g„(t,s)ds = e-*^°, (3.48) 

Jo 



where qa{t, s), s > 0, is the density of the stable process 5"(i), t > 0. 
We have that the probabilities 



P,^(t) = Pr{^^(5"(i)) = fcK^(0) = l} 

fc-l k 



(3.49) 



n^.E 



1 



j = i m=i lll=i,;^ 

fc-l fc 



1 m=l Ill=l,l^m. (-^i ~ ^™) "'" 

1 sir 



E^^i{~Xms'')qait, s)ds 



n^.E 



00 X 



j = l m=l l.\.l=ljyim\'^l ■^m) JO < > JO 

^^ A 1 sin^/TT r°^ r"-^ 

M 'i^iUlu^rni^i->^rn) ^ L r^'^ + 2r^ cosi^n + 1 

fc-l fc 



e ^ "^''qa{t,.s)dsdr 



e~*^'"''"dr 



Ra.E 



1 



1 ntl,;^r„ (^; - ^"0 "^ ^0 W" - + 2W = cos l^TT + 1 



j = l m 
fc-l fc 

n A, ^ — 

j = l m=l llz=l,;^m (A; — A,„) 



W 



^-1 



-iA^i^ -(1! 



dw 



1 



]Ee-*^^w„ 
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where Wa is a random variable with density 

/w„H = ^ ,^ , 7^ — , «;>0, 0<;.<1, (3.50) 

CCTT W^ Q + 2w Q COS UTT +1 

(first obtained by Lamperti ICI]). The density p.50|) coincides with the probability distribution of 

W. = (!)", (3.51) 

where S^, 82^ are independent stable random variables with Laplace transform 

Ee-^^^=e-^^ z>0,Q<v <l. (3.52) 

If a = 1, Ee^*"^"^" = Eyi{~\,-at'^) and p.49p are the state probabilities of a fractional pure birth process, 
while for a = v, Ee"*^™^" are the state probabilities of a pure birth process at time fWi, or, equivalently, a 
pure birth process at time t with rates XkWv- 

If we compare (|3.49p with (I2.55P , we can conclude that the process 

^''(5"(i))=^(T2,(5"(t))), (3.53) 

can be represented as 

^ (s'^'" (tWa)) , (3.54) 

\iQ < a <v <l. 

Remark 3.5. From formula p. 491) . when a takes the form, a — i^/2", n G N, we have 

(3.55) 

j = l m=l \.\.l=l,l=^m' ^' 



p^(i = / — If A, > —: dr. 



For n — ^ cx), we obtain that p'jl(t) — ^ 0, fc > 1, and pi (i) — >■ e *. T/iis shows that for n — >■ 00, either the 
population istantaneously explodes or does not produce offsprings with exponential probability. 

An alternative way of presenting the state probabilities p.49p is based on the Mellin-Barnes representation 
of the Mittag-Leffler function 

1 r+*°° r(2)r(i -z) ^ , , 

2m 7^_,^ r(/i - vz) 

with 1/ > 0, a; € C, |arg(-a;)| < tt (see Kilbas et al. I9|, page 44, formula (1.8.32)). 
In view of p.56p . we can write (|3.49p as follows 



"'" = n >■> E nt. „„.»,- A..) [^Jll ^P^(A-T'<i.,.(M).. (3.57) 



k—1 k ^ ^ /''7+'ioo 



tt\ \^ 1 1 r r(z)r(i -z)^ , f°" ,, , , , , 

II^.Et^I T^-^-.2m m-.z)^ - s'^^q^{t,s)dsdz, 

where in the last member of p.57p the Mellin transform of qa{t, s) appears. 
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The Mellin transform of the stable subordinator 5"(i), with Laplace transform 

reads 

l^/l-?7\ 1 



E(5"(i))''-i = / s'^-\^{t,s)ds ^ -r 
Jo a 



a y r(l -r;) 



By inserting p.59p into p.57p . we arrive at 

fe-i fe 



Pk 



w=n^.E 



1 



1 /'7+ioo 



j=l m=l l.\.l=ld^m\'^l ^ 



r(z)r(i-z)r(^z) 



1 ?ri=l lli = lJ7^m V^' ■^m) "'7-400 V ; V y 

1 



-1 tr2,l 



A„, t ' 



(0,1), (0,z/) 
(0,l),(0,i^/a),(0,z/) 



We examine now in detail the case a = i' in the next theorem. 
Theorem 3.5. We have the following distributions: 

1. Prl^Ms^it)) = k} = Pr{9i{T2,{s''{t))) - fc} = Pr{9i{tWu) = fc}, 

2. Pr{^ {s''{T2At))) = fc} = Pr{9i{tWi) = k], 



for fc > 1, t > 0, where 



W„ = 



SiY 



Sk 



and has distribution p. 501) . 
Proof. For fc > 1 we can write that 

PrK(r2,(5''(t))) = fc} 
fc-1 fc 



1 



/ e-^"-'Pr{T2,{s''{t))eds} 
Jo 

fT2^{z, s)q^{t, z)dzds 



'^mS 



j = l m=l Ill=l,l=£m y^l ^ 

fe-1 fe , 

j = l m=l [[l=lA=im i^l ^ ^^^) "^0 
K 1 A; -| /lOO 

Y[Xj Y^ —^ — - / E^.i{-\„,z'')q^{t,z)dz 

j = l m=l ll/=l,;#m l^i ^ ^m) "'O 



1 



fc-1 fc 

j = l m=l n/=l,;5^m i^l ~ ^m) '^ "'0 
fc-1 fc 



,00 /"OO 

smj/TT / / re 

dr 







r^^ I 0^1^ 



2r'' cos VTT + 1 



q^{t,z)dz 



n^^E 



1 



smz/^ /°° r''-ie-*('-^™'^)^ 



1 m=l ll/=l,;5^m (^' ^m) 

r , ,,. ^ sinvK r 

Pii-HXtr") = fc}. 



r 



/o TT r'^'^ + 2r'^ COS vn + 1 

= Pr{rAt(iWn = fc} = Pr{f^(iW,) = fc} . 



,21. _i_ 2r'' cos z/vr + 1 
dr 



dr 



(3.58) 



(3.59) 



(A,„i°) ""rfz (3.60) 



(3.61) 



(3.62) 
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This concludes the proof of the first resuh. In order to prove the second resuh we write 



1 



OO /'OO 



fc-l k 

j = l m=l Ill=l,l^m i^l ^ ^m) "'0 "'o 
fc— 1 k _. „oQ 

= n ^^ E —k r- / e-^'-^fT,^{z,t)dz 



(3.63) 



e '"^qv{z^s)fT2i,{z,t)dzds 



j = l m=l n/=l,;^m (^' ^ ^m) "'0 
byipU ^j^^ ^^ ^ _ ^ 



-Bi.,l(-A5^jt'') 



, fc-l fc 



1 



sini^TT r"^ ^e '^"*'' 



-dr 



J.J. J zi^ T-rfc (' \, _ \ ^ TT r^'' + 2r'^ cos t/TT + 1 

IJ j = l m=l 111=1, l=^m\'^l ■^m) 
Sill Z^TT 7* 

Pr{^(tr) = fc} 5^—— —dr = Pr{:At(iWi) = fc} . 

TT r^" + 2r'' cos vn + I 

n 

Remark 3.6. By slightly changing the above calculations, we arrive at the following result (compare with 

fe-i fc 



PrH (-S^m.W)) = fc} = n ^^- E TFT- 



1 



-Ee" 



-tA^lWi 



j = l m=l llz=l,;^m t'^' '^">) 

Remark 3.7. ^n alternative form of the distribution p.49p . /or a = i^, can 6e given as follows. 



Pr{^''{s''{t))=k} 

fc-i fc 



(3.64) 



(3.65) 



n^.E 



1 



1 



j = i m=i lll=i,;5^ 
fc-l fc 



1 1 



1 



1 



■=i ^^l^LM™(^'-^™)2*'^^ 



tA„ 



""''El (tA^e—^^) - e 



.-t^^-^dw 



tA^e" 



El (iAmC* 



ryXmt COS l^TT 



fc-l fc _. 

TT A,- y ^ i 



[e*^'""'"''''Ei(tA„ 



-27ri.'\ „iA^ sin y-n 



El (tA„e^-'^)] 



where the function Ei(z) — J ^-r-dt, \ argz < 7r|, is the exponential integral. 



3.4. Fractional pure birth process stopped at T2a{t) 

In this section we consider the process ^'^{T2a{t)), t > (see the discussion related to formula (I3.6p . As 
we did before, here we treat the case fc > 2. The state probabilities can be written as follows. 



P^ W = PrK''(T2.(t)) = fc [^"-(O) = 1} 

k 1 k ^ poo 

= n ^J- E TFfc TT TT / E,,,{-X^snPr{T2^{t) G ds} . 

j = l m=l [h=l,l^m {"^l ^ -^m) "'O 



(3.66) 
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The integral in ()3.66|1 can be further worked out by means of the Laplace transform: 

/ ^,.,i(-A™s'') / e-^'Pr {T2^{t) e ds} dt (3.67) 

Jo Jo 

= / ^.i(-A,ns")z"'^e-""°ds^z"-^ ^ ' = — . 

JO Z + A„i Z i- A,„ 

By taking the inverse Laplace transform of the above formula, we immediately obtain that 

/>oo 

£;..i(-A™s'')Pr{r2a(i)eds}= / ^,,i(-A™s'')i-"M^-_„,i-a(-i~"s)ds = S,„,i(-A™r"). (3.68) 

JO 

Therefore, the state probabilities for the process ^'^{T2a{t)), t > 0, result in the following form: 

k-l k 

pT^t) = n ^^ E TTk .. _. E,o.,i{-\ntn = PrW, k > 2. (3.69) 

j = l m=l lh=ld^rn (A/ A^) 

Note that the case fc = 1 can be treated in the same manner. We thus obtain the following equalities in 
distribution; 

^''(r2a(t))=rAt{r2,(r2«(i))}-^(T2,„(t)) = ^''"(t), t > 0. (3.70) 

Let now ?7„ = ni=i ^i^ where n € N, and Vi are n indices such that t'j € (0, 1] for I < i < n. Formula (|3.70p 
can be generalised as 

^{T2„,(T2„,(...r2,„(t)...))} = ^(r2„„(i)) = ^""(i), t>o (3.71) 

where 91"^" {t) is a nonlinear fractional birth process. 

Formula (j3.66p can also be worked out in an alternative way. In the following calculations we will make 
use of the integral representation 



vTit) = Pr{^''(T24t)) = k I ^''(0) = 1} (3.72) 

n ^^ E TTk 7T r- / ^-.1 (-A^s--) Pr{r2a(t) G ds} 

— 1 m=l ll;=i,/#m (A/ - Am) Jo 



3 

, fe-1 k 



Ha.E 



j=i m=i lli^ij^tm (A; A„ 
sinz/TT ^I'-i 



Jo TT r2'^ + 2r^ cos i/TT + 1 i ^"W ; 

M ' i'l nU^™ (A. - A„0 io - r-+2r^cos.. + l ^H 
Remark 3.8. iJy comparing formulae (j3.69p a77.rf (j3.72p . it is clear that the following expansion holds: 

sin i/TT f°° r""^ / IN 

E.o.Ai-Xmf'') = / ,, , ^ , — -Ba.i (-rA,^ni" ) dr, i^ G (0, 1], a G (0, 1]. (3.73) 

TT Jo ^ + ^r"-^ cos VTT + 1 \ / 



H^e give a direct proof of p.73p by applying the Laplace transform to both members. Of course 

fj.""^ + X 



OO ,,1/Q — 1 

e 



^*i?.a4(-A™t''")di = 1 , , ■ (3.74) 
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Then, we must calculate the twofold integral 



smi^TT I "~ ,,, /™ ^i' 1 



. 7o ^"Vo r^- + 2.-cos.. + l ^-^('-^^-^")-^^-^- ^'-''^ 



^ar 



I" J a T'^''' + 2?"'^ cos UTT + 1 



— iti f fi^ -\-rX 



/•oo /'oo a — 1 I/— 1 

-draw 



TT 7o 7o ^^'^ + 2r^ cos i/TT + 1 






Remark 3.9. A number of interesting relations follow from formula p.73p . 
T/ie following integral relation holds: 



smi^TT 



I' Jq r^ + 2r cos 



—E^,i (-r^A^r) dr (3.76) 

;os j/TT + 1 V / 



sin^/aTT /" 1 

e 



— ^""'dr. 



va Jq r^ + 2r cos i^an + 1 
A sort of commutativity is valid for p.73p ; 

Eo..A-^mtn = / 20^0 o TT^-i (-'■^™*i . ^ € (0, 1], a e (0, 1]. (3.77) 

TT Jo '' + 2?" cos aTT + 1 V / 

For a — 1 we recover, from p.73p . the integral representation of Mittag-Lejfler functions. By considering 
that 

^"^''^ ^ r^- + 2r-cosv^ + l' ^^'^^^ 

for V = 1, becomes a delta function with pole at r = \, we extract, from (I3.73p . an identity. 

Furthermore, it is worth noticing that formulae similar to (j3.73p can be derived by repeated applications 
of the same formula. For example we have: 

S.c./3,i (-A„,i''"'^) = ^ / ,^ , /^ -tEc^.i [~rX-^t^^] dr (3.79) 

^ ' TT Jo r'"^ + 2r'^ COS v-K + 1 \ / 

. n T^ — ?; r dw dr. 

Jo [r^" + 2r''cosvT: + l){w^°' +2w°'cosaT: + 1) 

Let Vi, 1 < i < n be n indices .such that for all I < i < n, vi G (0, 1], and let us denote rjn — Y[7=i ^i- ^^ 
general, for n > 2, we obtain that 

^^„,i(-A„,i''") (3.80) 

n ^i;--ri^ — -T ^-^1 --1-2^ -r- ■ ■ ■ r:ir a™-^ t- n ^-i 

^ r, ' + 2r^.' cos i/^ + 1 y \ / J-J[ 



1 1 



where ^^'Wi, 1 < j < n— 1, are independent random variables, each with distribution (|3.50p . wii/i a = 1 ond 
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